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NOMENCLATURE 
Coefficient  of  viscous  damping 

b 

Correlator  bandwidth 

D 

Dimensionless  damping  factor 

L 

Expected  value 

E 

Modulus  of  elasticity 

f 

Frequency  in  cycles  per  second 

u 

Undamped  natural  frequency 

Input  as  a function  of  space  and  time 

■f.Cn'b) 

Output  as  a function  of  space  and  time 

Transfer  function  or  generalized  response  function 

h(i-,l) 

I 

i 

Unit  impulse  response  function 
Logarithmic  decrement 
Imaginary  number,  \ - 

IX 


k 

Equivalent  spring  constant 

1 

Length  of  beam 

m 

Equivalent  mass  of  structure 

Generalized  mass 

Lateral  force  on  the  structure 

Position  vectors  on  the  surface  of  the  structure 

(ir,cls 

Elemental  surface  area 

Nonstationary  autocorrelation 

Autocorrelation  of  input  force  at  point  S 

Rio(r>^j^> 

) Cross-correlation  of  two  random  space-time  functions, 

and  positions  1"  and  S , 

re  spectively 

5(co) 

Power  spectral  density  as  a function  of  the  angular  fre- 
quency in  radians  per  second 

Time 

T 

Smoothing  time  or  time  length  of  the  sample  record 

X 


Lateral  displacement  of  the  structure 

Normal  mode  shape  of  motion  at  the  position  of  force 

YM) 

Normal  mode  shape  of  motion  at  the  accelerometer 
position 

oC 

Dimensionless  parameter  representing  the  autocorre- 
lation delay  X 

Dimensionless  parameter  representing  the  correlator 
smoothing  time  T* 

Generalized  damping  coefficient 

Time  variables 

Position  of  accelerometer 

K 

Natural  frequency  for  single  degree-of-f reedom  system 

^ tn 

Natural  frequency  of  the  m^  mode 

Damped  circular  frequency 

z 

CTkx 

Variance,  mean  square  variation  about  the  mean 

(K] 

V <T  ayji. 

Signal-to-noise  amplitude  ratio  of  the  correlator 

XI 


X 

Time  delay 

^..(r,u3) 

Power  spectral  density  of  the  input  force 

1 Cross-power  spectral  density  of  and  , 

re  spectively 

CO 

Angular  frequency  in  radians  per  second 

Natural  frequency  of  the  n^”  mode  in  radians  per  second 

-K- 

Denotes  complex  conjugate 

Xll 


CHAPTER  I 


INTRODUCTION 

One  of  the  criteria  for  the  mechanical  design  of  a structure 
involves  a consideration  of  the  structural  response  to  a discrete  fre- 
quency vibration  excitation.  In  the  case  of  a random  loading,  the  prob- 
lem becomes  more  difficult  and  time  consuming  since  we  now  must 
deal  with  a continuous  spectrum  of  vibration  excitation.  Before  it  is 
possible  to  predict  the  probable  response  of  structures  to  random 
loads,  one  must  first  define  the  loading  which  will  probably  occur  on 
the  structure. 

The  objective  of  this  study  is  to  experimentally  evaluate  the 
transfer  function  of  several  structures  under  stationary  environmental 
conditions.  To  attain  this  end,  the  random  pressure  from  a subsonic 
air  jet  has  been  regarded  as  a statistically  defined  load.  Once  the  load 
has  been  defined,  the  problem  reduces  to  the  determination  of  the 
transfer  function  and  then  to  the  prediction  of  the  probable  response. 

For  the  given  statistical  input,  it  is  possible  to  measure,  on 
a specimen,  the  transfer  function  as  a function  of  frequency  for  the 
points  of  interest.  From  these  measurements  on  the  structure,  it  will 
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then  be  possible  to  ascertain  the  validity  of  other  theoretical  repre- 
sentations of  the  same  structural  characteristics. 

A common  method  of  analysis  makes  use  of  several  normal 
modes  of  a coupled  spring-multiple  mass-dashpot  system  and  permits 
the  calculation  of  the  transfer  function  for  these  coupled  modes  of  vi- 
bration. The  larger  number  of  normal  modes  actually  present  in  the 
structure  may  not  allow  one  to  consider  such  a model  as  describing 
the  measured  curve  adequately. 


A single  degree-of-freedom  system  has  for  its  magnification 


factor 


H 


1 


(1.1) 


in  which 


■p  = the  impressed  frequency 


J = critical  damping  ratio 
k = equivalent  spring  constant 


hi  = mass 


The  response  to  an  excitation  at  a particular  frequency  is 

Response  = (Excitation  at  ) • (Transfer  function  at  -f^  ) 


(1.2) 
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where 


Er  = rms  excitation  at 

'a 

Hr  = transfer  function  at  -f 
W = rms  response 

To, 

For  several  discrete  frequency  excitations,  the  rms  response 


IS 


U = 


(1.3) 


and  may  be  extended  to  a continuous  spectrum  of  excitation  as 


u = 


XKhJJw 


(1.4) 


where  El^  is  the  power  spectral  density  of  the  input  force.  The  re- 
sponse of  a complex  structure  vibrating  in  all  its  modes  could  then  be 
evaluated  by  superposition  of  the  modal  responses. 

Aside  from  the  fact  that  the  results  for  several  modes  may 
not  give  an  accurate  description  of  the  response,  the  time  required  to 
go  through  the  above  process  is  not  justifiable.  With  this  in  mind,  the 
response  function  is  determined  for  a nonuniform  spectral  density  in- 
put. An  experimental  program,  based  on  a cross-correlation  theorem, 
is  undertaken  to  determine  the  transfer  function  for  several  structures 
and  in  one  case  is  compared  with  the  theoretical  modal  approach. 
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It  is  known  that  the  input-output  cross-correlation  of  a linear 
system  is  the  convolution  of  the  unit-impulse  response  and  the  input 
autocorrelation.  When  the  effective  bandwidth  of  the  input  spectral 
density  $(u)  is  not  white,  i.e.,  not  independent  of  frequency  over  a 
relatively  broad  frequency  interval,  it  is  not  an  easy  matter  to  obtain 
the  transfer  function.  The  experimental  determination  of  H(F,u) 
is  carried  out  by  working  exclusively  in  the  frequency  domain.  If  the 
input  force  power  spectral  density  had  been  white  noise,  the  cross- 
correlation of  the  output  acceleration  with  input  force  would  be  directly 
proportional  to  the  unit-impulse  response  [1].*  In  a block  diagram 
this  is  expressed  as  follows. 


A special  purpose  correlator  was  constructed  to  serve  this 
purpose.  The  method  of  cross-correlation  has  a distinct  advantage 
r other  conventional  methods  for  measuring  or  h("b) 


ove 


*The  numbers  in  brackets  refer  to  entries  in  the  List  of 
References . 
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Instead  of  applying  a step  input  and  recording  the  response,  one  can 
now  measure  the  transfer  function  under  stationary  conditions.  Thus, 
for  a stationary  random  input,  the  input  can  be  regarded  as  the  exci- 
tation for  the  unit-impulse  measurement.  This  is  valid  if  the  input 
noise  is  flat  and  wide  compared  to  the  bandwidth  of  the  system.  The 
procedure  is  to  measure  the  cross-correlation  R.  (KST)  and  R;;(sx) 


and  then  to  find 


by  a spacewise  F ourier  transformation  of 


and  However,  considerable  computation  is 


necessary  to  produce  results  of  acceptable  accuracy. 


CHAPTER  n 


APPLICATION  OF  CORRELATION  AND  SPECTRAL 
ANALYSIS  TO  THE  VIBRATION  OF 
ELASTIC  SYSTEMS 

A statistical  treatment  for  any  linear  structure  with  known 
normal  modes  excited  by  random  forces  of  given  power  spectral  den- 
sity of  excitation  has  been  carried  out.  The  mean  square  response 
can  be  evaluated,  and  the  probability  distribution  of  the  response  ex- 
ceeding an  assigned  value  is  then  available  under  the  assumption  of  a 
Gaussian  distribution  into  the  system.  The  response  of  these  systems 
has  always  been  determined  without  inclusion  of  spacewise  effects. 
They  have  been  studied  by  Van  Lear  and  Uhlenbeck  [2],  Eringen  [3] 
and  Thompson  [4]. 

The  basic  philosophy  of  the  treatment  of  random-in-time  phe- 
nomena has  been  given  by  Rice  [5].  The  restriction  that  the  noise  be 
a "stationary”  process,  and  furthermore  that  it  be  "Gaussian"  will  be 
assumed,  although  the  latter  assumption  is  not  needed  in  computing 
the  system  response  function.  The  fact  that  the  process  is  stationary 
Gaussian  allows  us  to  completely  specify  the  process  by  its  correlation 
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function  or  spectrum.  Furthermore,  an  ensemble  whose  sample  func- 
tions are  sums  of  sample  functions  from  Gaussian  random  processes 
is,  as  might  be  expected  from  the  central  limit  theorem,  the  ensemble 
of  a Gaussian  process.  Thus,  the  assumption  that  the  physical  dis- 
turbance of  the  structure  is  caused  by  innumerable  random  disturb- 
ances justifies  the  assumption  that  the  phenomena  be  Gaussian.  The 
assumption  of  a stationary  process  is  important  because,  otherwise, 
the  statistical  properties  may  change  in  time. 

Only  for  a linear  system  is  it  possible  to  relate  by  linear  op- 
eration certain  averages  of  the  output,  e.g.,  the  correlation  function, 
spectrum  and  average  power  to  the  correlation  function  or  spectrum  of 
the  input  process.  No  more  than  the  spectrum  of  the  input  need  be 
known  to  compute  the  corresponding  quantity  for  the  output  process. 
This  does  not,  by  any  means,  say  that  the  correlation  or  spectrum  de- 
scribes the  output  process  completely.  Unless  we  assume  that  the 
process  is  Gaussian,  the  spectrum  and  correlation  provide  little  in- 
sight in  determining  the  probability  distribution  of  the  response.  For 
a stationary  process  the  correlation  may  be  written  as 

^ (2.1) 

since  the  time  average  cannot  be  a function  of  time  and  can  depend 
only  on  time  differences.  Here,t  is  the  ensemble  average  and  R, 
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denotes  the  covariance  or  autocorrelation.  When  we  further  make  the 
assumption  that  the  process  is  ergodic,  we  can  also  say  that 


R(i,,0  - R,(-t)  = lim  i-  r'^x(t)x(i+i:)Jt  (2.2) 

In  the  above  we  have  interchanged  an  ensemble  average  for  a time 
average . 

It  should  be  mentioned  that  the  well-established  techniques  of 
spectral  analysis  developed  by  Rice  [5]  are  not  in  general  applicable 
to  nonlinear  systems,  although  the  analysis  can  be  extended  to  cover 
systems  containing  small  nonlinearity  [6],  A structure  with  nonlinear 
damping,  for  instance,  would  not  theoretically  permit  us  to  employ  the 
concepts  of  linear  harmonic  analysis  to  analyze  its  response  and  would 
consequently  require  the  determination  of  the  probability  density 
function. 

The  basic  relation  that  determines  the  response,  u , of  a 
beam  or  plate  is  the  following  differential  equation: 


m 


u-Htu4L(u.)  = 


(2.3) 


where 


L i 


s the  appropriate  differential  operator,  D is  the  coeffi- 


.t 


cient  of  viscous  damping  and  ^ denotes  the  forcing  function.  For  the 


plate  the  operator 


. L , is 
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L - C + 


(2.4) 


where  C is  a constant  depending  on  the  velocity  of  wave  propagation 
and  the  properties  of  the  plate. 

In  an  eigenvalue  approach  we  consider  solutions  of  the  Sturm- 
Liouville  equation 


L( 


m u + L u = o 


(2.5) 


satisfying  the  relevant  boundary  conditions.  The  normal  mode  solu- 
tions are  taken  in  the  form 


(2.6) 


where  is  the  natural  circular  frequency  of  the  n*"^  mode  whose 

shape  is  given  by  the  corresponding  X(p)  , The  solution  to  the 
original  differential  equation  can  be  assumed  to  take  the  form 


Substituting  (2.6)  into  (2.5)  gives 


- mu), 


^ L[v(f")'  = 

»-n  L t,  j 


(2.7) 


(2.8) 


The  solution  of  (2.3),  using  (2.7)  and  (2.8),  becomes 


(2.9) 
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We  next  use  the  orthogonality  property  of  normal  modes,  i.e., 


J V (r)  y^(r)  dp  = 0 , (2.10) 


and  for  m = Y1 


M=  Uv"  WJp- 

n t/  ^r\  ' 


(2. 11) 


where  denotes  the  generalized  mass,  and  assume  that  there  is  no 

damping  coupling  between  modes,  i.e.. 


JkV,(f)V_(F)JF.. 


n ' ' m 


(2.  12) 


and 


(2. 13) 


where  is  a generalized  damping  coefficient  (1/sec). 

On  multiplying  equation  (2.9)  by  and  integrating  over 

the  area,  this  equation  becomes 


(2. 14) 
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Taking  the  Fourier  transform  of  the  random  input  pressure  as 
and  the  output  as 

$f  / ii  W e''“^di  (2. 16] 

'n  ^ «» 

equation  (2.  14)  can  be  transformed  to 

where  is  the  natural  frequency  of  the  n^^  mode. 

To  introduce  the  power  spectral  density,  we  assume  that  the 
process  is  ergodic,  which  allows  us  to  write  the  spectral  density, S (u)), 
of  a single  function  from  the  ensemble,  i.e., 

W)  ^‘(tTu.)  ,2.18, 

the  asterisk  indicating  the  complex  conjugate.  Substitution  into  equa- 
tion (2.  17)  gives  the  output  power  spectrum  in  terms  of  the  input 


spectrum,  viz.. 
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where 


and 


0 = 


For  the  case  of  a rectangular  plate,  equation  (2.  19)  can  be  put  in 


another  form  [7], 


S.,( 


u 


rs.u) 

tl  tP 


V.(0  Vn(0 • 

/j’VJs)V,(5)S^(r,s,c)aA 
n '1.2, 

tn*t,i 


where 


is  a generalized  length  for  a beam  and  a generalized 


(2.  19) 


(2.20) 


(2.21) 


(2.22) 


area 


for  a plate,  and  the  represent  eigenvalues. 
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The  term  under  the  double  summation  sign  represents  the 
transfer  function  as  a function  of  space  and  frequency.  In  a subsequent 
chapter  this  term  is  specialized  for  the  case  of  a beam  by  deleting  one 
of  the  space  dimensions  and  then  summing  with  respect  to  five  modes. 


CHAPTER  m 


THEORETICAL  FORMULATION  OF  THE  PROBLEM 

The  purpose  of  this  study  was  to  attempt  to  apply  the  methods 
of  communication  theory  to  an  analysis  of  the  response  function  of 
structures  subjected  to  a loading  which  is  a function  of  space  as  well 
as  of  time.  The  derivation  of  the  cross-correlation  theorem  follows 
and  is  then  applied  to  the  determination  of  the  characteristics  of  a 
linear  system. 

In  a linear  system  subject  to  a random  force  input,  the  rela- 
tion between  the  input  force  autocorrelation  and  output  autocorrelation 
(of  either  displacement,  velocity  or  acceleration)  is 

f2^^(r,s,x)  = r KCr,s,t)<l»f  h(r,s/) R;i(r,s,T+v-<f  )ci\  d<T  (3.1) 

ioo  -oa 

where  ci^S  » dxj  and  the  integration  is  extended  over  the 

entire  region  of  the  structure.  Here  'if  and  (5  are  dummy  time 
variables. 

This  relation  is  found  by  using  the  definitions  of  input  and  out- 
put autocorrelation  functions  and  then  applying  the  convolution  integral. 
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(3.2) 


where  the  integration  is  extended  over  the  region  defined  by  the  input. 
Equation  (3.  1)  could  be  further  developed  in  order  to  obtain  the  unit 


sign,  its  evaluation  is  difficult.  The  unit  impulse  response,  or  Green's 
function,  is  more  readily  evaluated  by  considering  the  following  cross- 
correlation theorem  and  then  expressing  it  in  the  frequency  domain. 

For  a linear  structure  a theorem  exists  involving  the  input- 
output  cross-correlation,  the  input  autocorrelation,  and  the  unit- 
impulse  response,  h(r,x)  [8],  The  theorem  is  obtained  by  starting 
with  the  definition  of  the  input-output  cross-correlation 


impulse  response. 


, but  since  it  appears  under  the  integral 


(3.3) 


In  order  to  introduce 


into  this  expression,  we  make 


use  of  the  convolution  integral  for 


. Thus, 


(3.4) 
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where  scalar  force  field  defined  at  every  point  in  space. 

Inverting  the  order  of  integration,  we  have 


Since 


(3.6) 


equation  (3.5)  becomes 


(3.7) 


This  is  the  statement  that  relates  the  input-output  cross-correlation 
of  a linear  system  through  the  convolution  of  the  unit-impulse  response 
and  the  input  autocorrelation. 

To  express  this  in  the  frequency  domain,  we  take  the  Fourier 
transform  of  both  sides  of  the  equation,  viz., 

-itox  ,,  ioJX  I 

hCW)R..(nt-y)^^  0.8) 

-00  - 00  - 00 


17 


Introducing  the  Fourier  transform  pairs 


• 00 


H ( r,  s ,(0)  =J*  K (r,  s,x)  e'  ‘ <it 


4>nM=^/R,^^)e''“^ax 


and  letting 


J =x-■)^ 


^>io(''A“)  = ^ j'h(r,%,T)e '“^Jv/RaCnhe  JJ 


fi 


rA»)e‘"''a^^ 


- 00 


j*R|;W) 

.00 


>icoSi 

e dd 


(3.9) 


(3.  10) 


(3. 11) 


(3. 12) 


(3. 13) 


(3.  14) 


This  yields 
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$.Jr.s.u)  = H(r.s,u)  $,,(r.s,u) 


(3. 15) 


which  indicates  that  the  input-output  cross-power  density  spectrum  of 


the  input  power  density  spectrum. 

This  is  the  relation  used  in  computing  the  transfer  function. 
The  objective  here  is  to  indicate  how  this  basic  theorem  may  be  ex- 
tended to  measure  the  response  of  mechanical  systems  under  undeter- 
ministic  loads  that  are  functions  of  the  space  coordinates  as  well  as  of 
the  time.  This  method  also  intends  to  give  one  an  indication  of  the 
practicability  of  the  computations  involved  and  experimental  proce- 


a  linear  system  is  the  product  of  the  system  function 


and 


dures  required  to  carry  out  the  solution  of  the  problem. 


CHAPTER  IV 


EXPERIMENTAL  INVESTIGATION 

The  physical  models  studied  were  a cantilever  beam,  shown 
in  Figure  1,  and  a clamped  square  plate,  shown  in  Figure  2.  The  sup- 
port arrangements  are  depicted  in  these  figures.  The  beam  weighed 
90  grams  and  measured  7 inches  by  1 inch  by  5/16  inch.  It  consisted 
of  a lucite  plastic  inner  core  of  3/16  inch  thickness  cemented  to  an 
Armstrong  vibration  damping  material  1/8  inch  thick  on  either  side  of 
the  plastic.  The  damping  material,  identified  as  #881-658  FF,  de- 
creased the  time  necessary  to  obtain  acceptable  correlation  functions 
for  spectral  computation.  A similar  beam,  without  the  damping  ma- 
terial, called  for  time  delays  much  larger  than  the  afore-mentioned 
specimen.  The  5-inch  square  plate  weighed  approximately  210  grams 
and  was  constructed  of  the  same  material  as  the  beam.  The  total 
thickness  was  1/4  inch. 

The  experimental  flow  chart  is  shown  in  Figure  3.  The  input 
data  has  been  previously  recorded  on  a loop  spliced  from  a sample  of 
the  pressure  signal  sensed  by  a microphone  1/4  inch  in  diameter 
placed  in  a subsonic  air  jet  [9].  This  study  presented  statistical 
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Figure  1.  Beam  Test  Apparatus  and  Support 
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Figure  2.  Plate  Test  Apparatus  and  Support 


FM  Tape  recorder  no. 
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Figure  3.  Experimental  Flow  Chart  Showing  Arrangement  of  Equipment 
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representations  of  pressures  in  an  obstructed  field  [9].  It  was  this 
field  that  was  taken  to  represent  an  actual  service  load  on  the  struc- 
ture. The  transfer  functions  are  computed  by  cross-correlation  with 
this  force  under  stationary  conditions.  The  measurements  of  the 
pressure  were  made  with  two  crystal  microphones  1/4  inch  in  diame- 
ter, with  a frequency  range  of  20  to  90,  000  cycles  per  second.  The 
two  microphones  provided  simultaneous  measurement  of  pressure, 
at  a pair  of  points,  r and  S , on  the  surface  of  a rigid 
plate.  The  data  recording  system  is  shown  in  Figure  4. 


I Random  pressure 


Figure  4.  Recording  of  Random  Pressure 
from  a Subsonic  Air  Jet 
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A length  of  the  magnetic  tape  was  selected  in  which  the  proc- 
ess could  be  considered  to  be  acceptably  stationary,  and  then  the  tape 
was  spliced  into  a loop  which  corresponded  to  4 seconds  in  time  at  60 
inches/second.  This  time  is  called  the  recirculation  time  and  is  con- 
siderably greater  than  the  averaging  time  (time  constant)  of  the  corre- 
lator. In  Chapter  VI  the  effect  of  analyzing  a finite  length  loop  is  dis- 
cussed in  greater  detail.  The  length  of  loop  becomes  important  when  — • 

considering  the  spectrum  for  low  frequencies  because  the  resonances 
of  the  structure  on  which  the  vibration  is  measured  have  extremely 
narrow  bandwidths  at  this  end  of  the  spectrum. 

Typical  records  of  the  signal  and  its  autocorrelation  are 
shown  in  Figures  5,  6 and  7.  The  signals  were  recorded  simultane- 
ously on  a two-channel  instrumentation  tape  recorder. 

The  transfer  function  was  determined  for  one  point--in  the 
case  of  the  beam,  the  free  end,  and  for  the  plate,  the  center.  A co- 
ordinate system  was  established  so  that  position  0011  located  the  ac- 
celerometer. The  signal  was  detected  by  a Massa  accelerometer, 
model  M-142B.  This  unit  is  designed  so  that  it  can  be  mounted  on  the 
beam  through  a tapped  hole  in  the  base  of  the  unit.  The  response  of 
the  accelerometer  is  nearly  unity  because  of  the  high  frequencies  of 
crystals.  The  sensitivity  of  the  accelerometer  coming  out  of  the 
crystal  is  2.  88  mv/g.  The  additional  mass  of  the  accelerometer  is 
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Figure  5.  Oscilloscope  Record  of  Force 
and  Acceleration  of  Beam 


Figure  6.  Oscilloscope  Record  of  Force 
and  Acceleration  of  Plate 


Correlation 
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Figure  7.  Autocorrelation  of  Force  at  Position  13 
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assumed  to  be  negligible  in  comparison  to  that  of  the  beam.  It  is  rea- 
sonable to  suspect  that  the  accelerometer  mass  causes  greatest  error 
in  the  transfer  function  of  the  structure  at  the  position  0011.  This  is 
the  case  when  the  force  and  accelerometer  are  at  the  free  end. 

The  accelerometer  is  provided  with  an  18-inch  low-noise 
microdot  cable  that  is  attached  to  the  M-114B  preamplifier.  This  pre- 
amplifier is  designed  with  an  extremely  high  frequency  range  so  that  it 
may  be  used  in  conjunction  with  the  accelerometer  and  is  properly 
protected  and  isolated  from  mechanical  vibrations. 

After  amplification  of  the  output  of  tape  recorder  no.  1,  the 
signal  was  applied  to  an  attached  coil  (4  grams)  and  stationary  magnet 
driver  unit  as  shown  in  Figure  3.  The  coil-magnet  simulates  the  ran- 
dom force  at  the  designated  points.  Since  this  arrangement  acts  as  a 
shaker,  it  is  heavily  damped  by  the  power  amplifier  and  must  be  de- 
coupled to  prevent  the  coil-magnet  from  acting  as  a simple  dashpot. 
This  can  be  done  by  inserting  an  electrical  resistance  between  the 
power  amplifier  and  the  coil  as  shown  in  Figure  3.  The  driver  unit 
is  moved  alter  each  simultaneous  recording  of  the  force  and  the 
acceleration. 

In  the  study  of  the  beam's  dynamic  characteristics  the  mass 
of  the  coil-magnet  force-generator  is  neglected  as  was  done  for  the 
case  of  the  accelerometer  mass.  The  ratio  of  the  weight  of  the  coil 


28 


to  the  mass  of  the  beam  is  small  and  justifies  the  above  assumption. 
The  determination  of  the  autocorrelation  of  force  at  each  point  enables 
one  to  measure  the  disturbance  of  the  force  field  caused  by  the  move- 
ment of  the  beam.  This  is  witnessed  by  the  change  in  the  spectral 
density  of  the  force  (see  Figure  8).  At  first  it  was  thought  that  one 
could  compute  the  Fourier  transformation  of  the  autocorrelation  at 


This  cannot  be  done  if  one  wishes  to  obtain  accurate  results. 

The  driver  assembly  permitted  force -acceleration  correlation 
at  seven  positions,  spaced  1 inch  apart  along  the  beam.  (The  transfer 
functions  at  the  free  end,  corresponding  to  these  positions,  are  shown 
in  Figures  24  through  37,  Chapter  VI.  ) For  the  plate,  the  same  ar- 
rangement allowed  cross-correlation  at  nine  points,  spaced  approxi- 
mately 1 inch  apart  over  the  plate.  (The  transfer  functions  at  the  cen- 
ter of  the  plate  for  these  positions  are  shown  in  Figures  38  through  55, 
Chapter  VI.  ) 


does  not  appreciably  affect  the  dynamic  characteristics  of  the  beam, 
the  same  cannot  be  said  of  the  plate.  Because  of  the  many  resonances, 
plates  have  an  extremely  small  driving-point  impedance.  Accurate 
measurements  are  difficult  to  obtain,  particularly  at  the  higher  fre- 


any one  point  and  subsequently  use  this  in  the  computation  of 


Although  the  additional  mass  of  shaker  and  accelerometer 


quencies  where  the  mass  impedance  of  the  coil  exceeds  the  driving- 
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Figure  8,  Structural  Feedback  Effect  on  Spectral  Density  Input 
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point  impedance  of  the  plate.  As  in  the  case  of  the  beam,  a resistance 
was  placed  in  series  with  the  coil-magnet  arrangement  in  order  to  re- 
duce the  feedback  from  the  coil  to  the  amplifier  as  shown  in  Figure  3. 


Correlation  Measurement 

The  exact  determination  of  the  autocorrelation  function  or 
power  spectral  density  would  require  an  infinitely  long  sample.  In 
addition,  the  sample  must  be  assumed  to  be  stationary.  A single  time 
function  of  the  acceleration  data  was  cut  into  loops  of  various  smooth- 
ing times,  T , and  the  correlation  functions  corresponding  to  these 
loops  compared.  For  each  loop  the  correlation  was  plotted  and  the 
Fourier  transformation  evaluated.  Since  this  showed  that  there  was 
no  significant  change  in  the  spectral  density,  the  samples  could  then 
be  considered  as  approximately  stationary  for  the  short  intervals  of 
time.  A signal  whose  averages,  when  performed  on  successive  sam- 
ples, show  no  change  with  time  is  called  a stationary  signal.  Thus  the 
correlation  can  be  written  as 


(4.  1) 


Since  the  smoothing  time  is  fixed,  the  correlation  curves  are  no  more 
than  estimates  of  the  true  correlation.  The  normalized  value  of 
is  obtained  by  dividing  the  time  average  of  the  two  signals. 
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and  , by  the  rms  value  of  each  of  these  signals.  The  normal- 

ized correlation  takes  the  following  form  for  electronic  computation: 


Rio(nt)=  (4.2) 

The  length  of  record,  T , was  chosen  as  4 seconds.  The  lag,  X , was 
taken  at  values  less  than  0.5  per  cent  of  the  length  of  the  sample. 
Theoretically,  the  computation  of  RJnx)  is  valid  only  for 
where  is  the  maximum  desirable  lag.  The  block  diagram  for 

this  operation  is  displayed  in  Figure  9. 


DC  voltmeter 


Figure  9.  Flow  Diagram  for  the  Evaluation  of  Correlation 
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The  acceleration,  , and  force,  f.(ni)  . 

are  played  back  simultaneously  from  channels  1 and  5,  respectively, 
of  the  Ampex  FR-1100  Recorder /Reproducer  system.  The  recorder 
is  equipped  with  seven  channels.  The  appropriate  time  shift  was  in- 
troduced by  an  audio  delay  line  with  a time  delay  range  of  0 to  2000 
microseconds  on  either  of  the  two  signals.  The  two  signals  were 
multiplied  with  a Philbrick  Multiplying/Dividing  component  with  a 
15  kc  bandpass.  The  time  average  of  the  multiplier  output  is  formed 
by  a dc  vacuum-tube  voltmeter.  The  rms  values  of  the  respective 
signals  were  found  by  the  same  procedure  as  above,  except  that  the 
time  shift,  X , was  not  introduced. 

It  may  be  observed  from  Figure  10  that  the  time  delay,  X , 
is  extended  to  10  milliseconds.  This  large  delay  is  made  possible  by 
delaying  a signal  for  2000  microseconds  and  then  recording  it  on 
another  channel.  This  can  be  repeated  as  often  as  desired,  since  the 
fm  track  will  erase  any  previous  signal  as  it  records.  Thus,  a 2000- 
microsecond  delay  in  the  force  is  recorded  on  channel  5,  and  at  the 
same  time  the  acceleration  is  recorded  on  channel  7 without  delay. 

In  the  computation  of  the  correlation  function,  the  interchannel  time 
displacement  error  must  be  talcen  into  account.  This  is  caused  by  the 
misalignment  of  the  track  heads.  Misalignment  produced  a time  delay 
error  as  large  as  25  microseconds. 


Correlation 
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A comparison  of  correlation  functions  plotted  for  different 
smoothing  times  showed  that  the  quality  of  the  spectral  density  esti- 
mate does  not  increase  insofar  as  the  computation  of  the  transfer 
function  is  concerned. 

The  acceleration  and  force  were  recorded  for  20  seconds 
at  a tape  speed  of  30  ips.  These  records  were  cut  and  spliced  into 
4- second  loops.  During  the  playback  process,  the  tape  speed  was 
increased  to  60  ips. 


Power  Spectral  Measurement 

It  was  previously  mentioned  that  a theoretically  accurate 
determination  of  the  power  spectral  density  requires  an  analysis  of  a 
stationary  random  phenomenon  during  an  infinite  time  interval.  Since 
this  cannot  be  done,  the  results  obtained  are  estimates  of  the  power 
spectrum,  even  though  the  operations  involved  in  the  process  are  exact. 

After  the  correlation  function  was  obtained,  the  Fourier 
transformation  required  for  establishment  of  the  power  spectral  den- 
sity was  performed  on  the  IBM  650  computer.  The  computation  is 
incorporated  in  the  following  formula: 


■/ 


R (t)»R 

mn 


0 •- 


nm 


lx) 


cosiiT-txdxt 


•'/ 


O 


R.  W-R  Wlsinin^xch:  (4.3) 

nm  mn 
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The  integral  is  divided  into  a finite  series  over  an  interval,  AT  , 
where 

= k At  k=  0^1, L,.. . . 


The  final  numerical  computation  for  the  cross-power  spectral  density 
become  s 


= £.-C 


+iAt^(R^^-R^Psm(o-C| 


(4.4) 


kri 


v-6 


where  At  is  taken  as  40x10”'^  seconds.  When  X*y  , $ XX 

comes  the  power  spectral  density, 


(05) --2  AX 


% R 


^^Cos 


(4.5) 


A typical  set  of  curves  computed  from  the  above  formulas  for 
several  positions  on  the  beam  is  shown  in  Figures  11  through  15.  Simi- 
lar curves  are  indicated  in  Figures  16,  17  and  18  for  the  plate. 

It  should  be  mentioned  that  the  cross-power  spectrum  cannot 
be  interpreted  as  a transfer  of  power  between  different  parts  of  the 
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Spectral  density,  lO""^  sec 


Figure  11,  Cross-Power  Spectral  Density  of  Force-Acceleration 
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Spectral  density,  10"'^  sec 


Figure  12.  Cross-Power  Spectral  Density  of  Force-Acceleration 
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13.  Cross-Power  Spectral  Density  of  Force-Acceleration 
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Figure  14.  Cross-Power  Spectral  Density  of  Force-Acceleration 
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Power  spectral  density,  lO""^  sec 


Figure  15.  Power  Spectral  Density  of  Input  Force 


Correlation 
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Spectral  density,  10"'^  sec 


Figure  18,  Cross-Power  Spectral  Density  of  Force-Acceleration 
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structure  or  the  spectrum  of  power  taken  by  the  structure.  Although 
the  cross-correlation,  R jty  , after  transformation,  has  the  same 
dimensions  as  power,  it  does  not  yield  a spectrum  of  power.  The  in- 
put force  and  output  acceleration  are  cross-correlated  at  different 
locations  on  the  structure  merely  for  transfer  function  determinations. 


CHAPTER  V 


NORMAL  MODE  DETERMINATION  OF  THE 
TRANSFER  FUNCTION 


An  attempt  has  been  made  to  compare  the  beam  transfer  func- 
tion as  determined  from  the  cross-correlation  theorem  with  a model 
that  is  established  in  terms  of  a simple  harmonic  motion. 

Referring  to  the  analysis  on  page  12,  equation  (2.22),  the  re- 
sponse $^(r,S^Co)  can  be  written  as 


$i  (^.  <jr^ds. 


This  can  be  rewritten  as 


(5.  1) 


'it 


1^ 


(5.2) 
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Equation  (5.Z)  is  put  in  the  following  form,  where  the  conjugate  is  de- 
noted by  an  asterisk  (*): 


The  complex  frequency  response  is 


m 


V™WV„W  [ 

.ft 

(5.4) 


where  D is  a dimensionless  damping  factor. 

In  this  form  the  transfer  function  is  readily  seen 

to  be  a superposition  of  modes,  each  considered  to  be  a single  degree- 
of -freedom  system  with  a particular  natural  frequency  and  damping 
characteristic  as  shown  in  Figures  19  and  20,  Equation  (5.4)  ex- 
presses the  transfer  function  for  displacement  of  the  beam.  To  find 
the  transfer  function  for  acceleration,  we  can  make  use  of  the  follow- 
ing expressions  for  derivatives  of  correlation  functions.  Starting  with 
the  definition  of  the  cross-correlation. 


(5.5) 


and  differentiating  under  the  integral  sign,  we  have 
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1 - (f/fn)^ 


[1  - (f/fn)^]  + 4D^(f/fn)^ 


Figure  19.  Real  Part  of  Beam  Response  Represented  as  a 
Single  Degree-of-Freedom  System 
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Figure  20.  Imaginary  Part  of  Beam  Response  Represented  as  a 
Single  Degree-of-F reedom  System 
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+ T 


-T 


(5.6) 


where  ^ Ko^)  “^o 


. A second 


differentiation  under  the  integral  sign  gives 


T 


-T 


Thus  the  second  derivative  of  the  cross-correlation  of  the  force  and 
displacement  is  another  cross-correlation  between  force  and 
acceleration. 

To  find  the  corresponding  relation  in  the  frequency  domain, 

we  write 


(5.8) 


and  differentiate  both  sides  twice;  thus, 


(5.9) 


Taking  the  inverse  transform,  we  obtain 
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2TT 


(5. 10) 


Since 


H(n“) 


1,0 

u 


the  corresponding 


H(r,co) 


for  acceleration  is  given  by 


(XCC. 


(5. 11) 


(5. 12) 


The  theoretical  and  experimental  transfer  functions  are  plotted  for  the 


acceleration  response  (see  Figures  24  through  56,  Chapter  VI). 

TABLE  1 


EXPERIMENTALLY  DETERMINED  BEAM  PARAMETERS 


Mode 

Tu 

Measured 
F requency 

Ratio  to  First 
Mode  Frequency 

Damping 
Factor  (D) 

1 

24.2 

1.  0 

0. 0442 

6.64 

2 

158.  0 

6.  5 

0. 0446 

44.  10 

3 

430.  0 

17.  9 

0.  0450 

109. 50 

4 

910.  0 

37.  9 

0.  0500 

285.00 

5 

1450.  0 

60.  4 

0.  0519 

473.  00 

6 

2250. 0 

93.  9 

0. 0555 

734. 00 
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Table  1 lists  the  natural  frequencies  found  experimentally. 
A Tektronix  dual-beam  oscilloscope  was  used  in  conjunction  with  a 


input  and  output  acceleration.  In  Table  1 the  first  six  modes  are 
listed. 

To  determine  the  damping  of  each  mode,  use  was  made  of  a 
Dumont  oscillograph  record  camera.  The  beam  is  excited  at  its  natu- 
ral frequency  with  a steady- state  sinusoidal  input.  At  a predetermined 
time,  a circuit  was  opened  which  triggered  the  sweep  on  the  oscillo- 
scope, and  this  permitted  a record  of  the  trainsient  vibration  to  be  ob- 
tained. Figures  21  and  22  show  a record  of  a transient  at  a frequency 
of  62  and  940  cycles  per  second,  respectively.  To  determine  the 
damping  coefficient  for  a particular  mode,  it  is  assumed  that  the  beam 
acts  as  a single  degree-of-freedom  system.  This  is  approximately 
true  for  a natural  vibration.  Thus,  the  differential  equation  that  fol- 
lows defines  the  model: 


wide-range  oscillator  to  produce  the  Lissajous'  figures  between  the 


(5. 13) 


where 


lateral  displacement 
K.  = natural  frequency 


1^.  - 


= damping  factor. 


•V 
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Figure  21.  Oscilloscope  Record  of  Beam  Transient 
for  Second  Mode 


Figure  22.  Oscilloscope  Record  of  Beam  Transient 

for  Third  Mode 
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The  solution  of  this  equation  is  written  as 


-6  ^ 

cos(^i^-cc) 


(5. 14) 


If 


we  denote  by  ^ ^ ..^c^^the  successive  maximum  displacements 


on  the  same  side  of  the  axis,  the  ratio  of  the  amplitude  to  the  ^ 

amplitude  will  be  given  as 


(5. 15) 


where 


and 


T ^ 


Here  — K 41-D*"  represents  the  circular  frequency  in  the  presence 
of  damping. 

The  logarithmic  decrement  I is  found  from 


I- 


hn 


^k  _ 


k+i 


ZTT  D 

h'D’’' 


(5.16) 


Figure  23  indicates  the  determination  of  V for  the  second  mode  at  a 
frequency  of  158  cycles  per  second.  Values  of  successive  amplitudes 
are  read  directly  from  the  oscillograph  records,  Figures  21  and  22. 
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Figure  23.  Plot  of  Logarithmic  Decrement  for  the 
Second  Mode  of  Vibration 
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V 


The  slope  of  the  curve  gives  the  damping  coefficient  . The  value 


not  assume  that  the  damping  coefficient  is  constant  for  all  modes,  as 
is  usually  the  case.  Table  1 lists  the  damping  coefficients  for  each 
respective  mode. 

The  characteristic  functions  representing  the  normal  modes 
of  a cantilever  beam  have  been  computed  by  Young  and  Felgar  [ 10] . 
Tables  2,  3,  4,  5 and  6 show  a sample  computation  of  at 

frequencies  of  500,  1000  and  1500  cycles  per  second.  In  Chapter  VI 
the  results  of  the  modal  solution  are  represented  by  the  dotted  curve. 


of  is  determine 


d by  multiplying  by  the  natural  frequency  in 


radians  per  second.  The  following  calculation  of 


doe  s 


J) 


f \ -1>'^ 


REAL  (R)  AND  IMAGINARY  (1)  PARTS  OF  THE  TRANSFER  FUNCTION 
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CHARACTERISTIC  MODE  SHAPES  FOR  CANTILEVER  BEAM 
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TABLE  5 

CONTRIBUTION  OF  THE  N*^^  MODE  TO  THE 
TRANSFER  FUNCTION 


Mode 

f/f„ 

Rn 

In 

f 

= 500 

1 

10.416 

-0. 00930 

+ 0.  000078 

2 

1.5822 

-0. 6594 

+ 0.  0624 

3 

0.5814 

+ 1. 5010 

+ 0.  1187 

4 

0. 2747 

+1. 0807 

+ 0.  0321 

5 

0. 1724 

+ 1. 0303 

+ 0.  0190 

f 

= 1000 

1 

20.833 

-0. 00231 

0 

2 

3. 1646 

-0. 1109 

+ 0.  00350 

3 

1. 1628 

-2. 6095 

+ 0.  7755 

4 

0.5494 

+ 1.4235 

+ 0.  11202 

5 

0. 3448 

+ 1.  1331 

+ 0.  04611 

f 

= 1500 

1 

31. 25 

-0. 0010 

0 

2 

4.7468 

-0. 0464 

+0. 00092 

3 

1.7442 

-0.48678 

+ 0.  037417 

4 

0. 8242 

+2. 9250 

+ 0.  75170 

5 

0.  5172 

+ 1.  3578 

+ 0.  09970 
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TABLE  6 

REAL  AND  IMAGINARY  COMPONENTS  OF  THE  TRANSFER 
FUNCTION  FOR  ACCELERATION  AND  DISPLACEMENT 
CONSIDERING  PARTICIPATION  OF  FIVE  MODES 

OF  VIBRATION 


Position 

HxlO"^ 

f = 500 

f = 1000 

f = 1500 

0011 

Displacement 

Real 

Imaginary 

-22. 555 
+ 1.  315 

-7.753 
+ 0.  875 

-2. 5526 
+ 0.  100 

0011 

Acceleration 

Real 

Imaginary 

+5. 6388 
-0. 3288 

+7.753 

-0.875 

+5. 743 
-0. 225 

0012 

Displacement 

Real 

Imaginary 

-12. 220 
-0.  030 

+ 0.  041 
-0. 858 

-1.  102 
-0.  173 

0012 

Acceleration 

Real 

Imaginary 

+ 3.  055 
+ 0.  0075 

-0. 041 
+ 0.  858 

+2. 480 
+ 0.  389 

0013 

Displacement 

Real 

Imaginary 

-1.  157 
-0.  906 

+ 3.  265 
-1.  153 

+ 0.  494 
-0.  063 

0013 

Acceleration 

Real 

Imaginary 

+ 0.  2893 
+ 0.  2265 

-3. 265 
+ 1.  153 

-1.  112 
+ 0.  142 

0017 

Displacement 

Real 

Imaginary 

+ 1.  148 
-0.  047 

-0.771 
+ 0.  227 

-0. 279 
-0.  236 

0017 

Acceleration 

Real 

Imaginary 

-0.  2870 
+ 0.  0118 

+ 0.  771 
-0. 227 

+ 0.  628 
+ 0.  081 

CHAPTER  VI 


PRESENTATION  AND  DISCUSSION  OF 
EXPERIMENTAL  RESULTS 

It  is  evident  from  this  study  that  precise  data  reduction  of  a 
random  vibration  with  an  irregularly  shaped  spectrum  of  force  is  a 
slow  process,  although  it  can  be  made  automatic.  The  time  needed  to 
produce  a correlation  curve  is  approximately  2 hours.  The  raw  vi- 
bration data  are  obtained  on  magnetic  tape,  and  in  preparing  for  a 
harmonic  analysis  a brief  portion  of  the  tape  is  spliced  into  a loop. 

The  loop,  when  played  back,  has  a definite  recirculation  frequency 
even  though  the  actual  measurement  is  assumed  random.  In  order  to 
estimate  the  error  imposed  by  the  finite  integration  time,  we  may  as- 
sume that  the  correlation  curves  shown  in  Figures  10  and  17  are  ex- 
ponential-cosine functions  and  are  obtained  from  a normal  random 
process.  It  can  be  shown  that  the  mean  square  variation,  V , 
about  the  mean  value  of  a set  of  measurements  is  given  by  the  follow- 
ing formula  [ 1 ] : 

where  k denotes  the  different  ensemble  members. 
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With  the  assumption  of  an  exponential  correlation,  this  re- 


duces to 


where  the  dimensionless  parameters  oL  and  are 

o(  = irBt 
ttBT 

This  is  an  exact  expression  for  all  oC  and  ^ for  which  the  delay  X 
is  less  than  the  smoothing  time  T*  . For  the  loop  under  consideration, 
the  smoothing  time  is  approximately  4 seconds,  and  the  bandwidth  of 
the  correlator  is  about  15,000  cycles  per  second,  considered  as  the 
bandpass  of  the  multiplier.  Thus, 

TT  BT  = TT(l5,0O0)(4)i  2 Xlo®  (6.3) 

This  value  of  y is  sufficiently  large  [ 1,  p.  269]  to  justify  use  of  the 
formula  developed  by  Bendat  [ 1,  p.  268],  i.  e.  , 


XX 


i<k<4 


(6.4) 


The  expected  rms  error  in  a set  of  measurements  of  the 
autocorrelation  function  at  any  point  along  the  autocorrelation  curve 
for  is  =0-004-5  [1,  Table  7.2,  p.269]. 
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For  all  values  of  T and  T.  , the  output  signal-to -noise  amplitude 
ratio  of  the  correlator  is  defined  as  the  ratio  of  the  expected  mean 
value  to  the  standard  deviation,  that  is, 

f_R_]  = (6.5) 

\ /o.V,  ) 

For  an  exponential  correlation  curve,  this  becomes 


/ R ] — 

e 

V ^ ioui. 

(6.  6) 


Since  '9  = ^ X iO  and  is  thus  considered  large,  we  can  substitute 
(6.4)  into  (6.6)  for  . Thus 


(6.7) 


for  large  y . 

This  formula  gives  a signal-to-noise  ratio  of  1.5  for  cC  = 5. 
It  shows  that  the  accuracy  in  the  correlation  measurement  drops  off 
quickly  for  values  oL  greater  than  5.  This  result  indicates  that  the 
output  signal-to-noise  ratio  will  be  less  than  1.5  for  a time  delay 
greater  than  100  microseconds.  The  exponential  correlation  function 
represents  an  upper  bound  in  evaluating  the  statistical  error 
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measurement  of  the  correlation  function.  It  is  reasonably  safe  to  say 
that  the  errors  for  an  exponential-cosine  correlation  will  be  much  less 
than  those  obtained  above  and  will  not  give  erroneous  results  for 

In  determining  the  autocorrelation  at  different  locations  along 
the  structure,  it  was  found  that  the  system  disturbs  the  force  field. 

This  was  shown  in  Figure  8 by  the  comparison  of  two  spectral  density 
curves  of  input  force.  The  autocorrelation  must  be  evaluated  at  every 
point  in  order  to  account  for  the  disturbance  created  by  the  structure. 

The  resulting  transfer  functions  for  a cantilever  beam  and 
clamped  plate  are  shown  in  Figures  24  through  56.  The  graphs  show 
the  amplitude  and  phase  as  a function  of  frequency  in  cycles  per  sec- 
ond. Each  graph  indicates  the  spacewise  contribution  to  the  transfer 
function.  In  Figures  24  through  30  the  dotted  curves  represent  the 
analysis  of  Chapter  V for  the  modal  vibration.  At  position  0017  the 
theoretical  solution  deviates  considerably  from  the  experimental  curve. 


fixed-free  beam.  Because  the  experimental  beam  is,  of  course,  not 
completely  fixed,  it  is  to  be  expected  that  the  in  the  proximity  of 

the  support  do  not  give  an  accurate  description  of  the  experimental 
characteristic  shapes.  For  points  far  from  the  fixed  end  the  results 
agree  favorably  with  an  assumed  fixed-free  beam.  The  partially 
clamped  end  has  little  effect  near  the  latter  region. 


The  normal  modes 


introduced  into  the  solution  are  for  a 
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Magnification  factor 


F requency,  cps 


Figure  24.  Modulus  of  Acceleration  Transfer  Function  at  0011 
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Magnification  factor 


Figure  25.  Modulus  of  Acceleration  Transfer  Function  at  0012 
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Magnification  factor 


Figure  26.  Modulus  of  Acceleration  Transfer  Function  at  0013 
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Magnification  factor 


Figure  27.  Modulus  of  Acceleration  Transfer  Function  at  0014 
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Magnification  factor 


Figure  28.  Modulus  of  Acceleration  Transfer  Function  at  0015 
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Figure  29.  Modulus  of  Acceleration  Transfer  Function  at  0016 
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Magnification  factor 


Figure  30.  Modulus  of  Acceleration  Transfer  Function  at  0017 
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Figure  31,  Phase  Angle  of  Acceleration  Transfe 


r Function  at  0011 
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Figure  32.  Phase  Angle  of  Acceleration  Transfer  Function  at  0012 
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Figure  33.  Phase  Angle  of  Acceleration  Transfer  Function  at  0013 
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Figure  34.  Phase  Angle  of  Acceleration  Transfer  Function  at  0014 
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Figure  35.  Phase  Angle  of  Acceleration  Transfer  Function  at  0015 
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Figure  36.  Phase  Angle  of  Acceleration  Transfer  Function  at  0016 
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Figure  37,  Phase  Angle  of  Acceleration  Transfer  Function  at  0017 
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Figure  38.  Plan  View  of  Plate 
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Figure  39.  Modulus  of  Acceleration  Transfer  Function  at  0011 
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Figure  40.  Modulus  of  Acceleration  Transfer  Function  at  0012 
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Figure  41.  Modulus  of  Acceleration  Transfer  Function  at  0013 
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Figure  42.  Modulus  of  Acceleration  Transfer  Function  at  0021 
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Figure  43.  Modulus  of  Acceleration  Transfer  Function  at  0022 
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Magnification  factor 


Figure  44.  Modulus  of  Acceleration  Transfer  Function  at  0023 
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Figure  45.  Modulus  of  Acceleration  Transfer  Function  at  0031 
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Magnification  factor 


Figure  46.  Modulus  of  Acceleration  Transfer  Function  at  0032 
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Magnification  factor 


Figure  47.  Modulus  of  Acceleration  Transfer  Function  at  0033 
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Figure  48.  Phase  Angle  of  Acceleration  Transfer  Function  at  0011 


Phase  angle,  rad/sec 


Figure  49.  Phase  Angle  of  Acceleration  Transfer  Function  at  0012 
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Figure  50.  Phase  Angle  of  Acceleration  Transfer  Function  at  0013 


95 


Figure  51.  Phase  Angle  of  Acceleration  Transfer  Function  at  0021 
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Figure  5 2.  Phase  Angle  of  Acceleration  Transfer  Function  at  0022 
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Figure  53.  Phase  Angle  of  Acceleration  Transfer  Function  at  0023 
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Figure  54.  Phase  Angle  of  Acceleration  Transfer  Function  at  0031 
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Phase  angle,  rad/sec 


Figure  55.  Phase  Angle  of  Acceleration  Transfer  Function  at  0032 
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Phase  angle,  rad/sec 


Figure  56.  Phase  Angle  of  Acceleration  Transfer  Function  at  0033 
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An  experimental  study  of  plates  for  determination  of  their 
transfer  function  is  much  more  difficult  than  would  be  expected.  Be- 
cause of  the  large  number  of  modes  in  which  it  can  vibrate,  the  plate 
exhibits  a very  small  driving  point  impedance.  Consequently,  the 
mass  impedance  of  the  coil  and  accelerometer  have  an  appreciable 
effect  on  the  measurements.  The  transfer  function  for  the  higher  fre- 
quencies is  therefore  questionable.  In  evaluating  the  transfer  function 
for  points  other  than  the  center,  it  will  be  necessary  to  drive  the  plate 
over  its  entire  surface.  The  same  techniques  as  before  are  involved, 
but  the  labor  is  greatly  increased.  For  the  center  of  the  plate,  the 
change  in  the  transfer  function  as  a function  of  space  is  most  easily 
seen  through  variation  of  its  phase.  With  increasing  distance  from 
the  center  of  the  plate,  the  space-phase  deviations  between  the  various 
modes  become  greater.  The  magnitude  of  the  transfer  function  re- 
mains relatively  constant  with  distance  from  the  center.  Since  the 
damping  is  quite  small,  the  system  is  likely  to  vibrate  with  the  same 
amplitude  all  over  its  area.  If  the  damping  were  larger,  the  space- 
amplitude  part  of  the  transfer  function  would  decrease  very  rapidly 
with  the  distance  from  the  center  of  the  plate. 

Studies  of  random  vibrations  of  continuous  structures  have 
been  carried  out  for  distributed  stationary  excitation,  uncorrelated  in 
both  space  and  time.  Thomson  and  Barton  [4]  indicate  the  solution 
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for  beams  when  the  excitation  is  confined  to  a single  point.  The  pres- 
ent study  allows  one  to  consider  spacewise  correlated  excitation. 

Once  the  transfer  function  is  evaluated  for  the  desired  points  of  inter- 
est, the  response  to  any  correlated  loading  can  be  determined.  The 
transfer  functions  remain  valid  for  the  case  of  a nonhomogeneous  load 
field  and  can  be  used  to  evaluate  the  response  as  long  as  the  system  is 
linear.  The  cross-correlation  technique  for  measuring  has 

the  advantage  that  one  cam  make  measurements  when  the  system  is  in 
normal  operation  under  service  loads.  This  is  clearly  am  advantage 
since  one  does  not  have  to  excite  the  structure  with  special  types  of 
forces  to  obtain  the  structure's  transfer  fvmction,  as  is  done  in  the 
case  of  conventional  sinusoidal  testing. 


103 


Figure  57.  Complete  Assembly  of  Experimental  Equipment 
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